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We generalize the Kohn's theorem to the case of parabolic three-dimensional quantum dots in 
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magnetic fields of arbitrary direction. We show numerically that the exact resonance frequencies in 
the magneto-optical absorption of these dots are reproduced by the adiabatic time-dependent local 
spin density approximation theory (TDLSDA). We use TDLSDA to predict spin density excitations 
in the dots. 

PACS numbers: 71.15.Mb, 73.21.La, 73.22.Lp, 85.70.Sq 
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I. INTRODUCTION 



Longitudinal dipole excitations in quantum dots (QD's) under applied magnetic fields 



have been observed in photo-absorption experiments in the far infra-red region, and 
more recently in Raman scattering experiments. ^, [4] These modes are excited by the dipole 
operator D = fj, and the resonance frequencies in the magneto-optical absorption spec- 
trum of a QD with asymmetric parabolic confinement potential are found to be independent 
of the electron-electron interaction and of the number of electrons in the dot, and they coin- 
cide with the single-electron transition frequencies. This statement, known as the generalized 
Kohn theorem, was demonstrated some years ago by Peeters|5| for a two-dimensional QD 
under an external magnetic field (B) perpendicular to the plane of motion of the electrons. 
Peeters gave the following explicit formula for the resonance frequencies: 



<2 , 



where u c = eB/mc is the cyclotron frequency, and uj x , u y are the confinement frequencies 
in the x and y directions. Time-dependent theories as time-dependent Hartree-Fock, and 
adiabatic time-dependent local spin density approximation (TDLSDA) fulfil the generalized 
Kohn theorem, p, Q, ^ and thus any numerical implementation of these methods should 
reproduce the exact result of Eq. (HI). Viceversa, static mean field theories like Hartree, 

Q 

Hartree-Fock and Kohn-Sham theories violate the theorem. [9( 

In this paper we generalize the result given in Eq. $1} to the case of a parabolic three- 
dimensional QD in a magnetic field of arbitrary direction, and present and discuss adiabatic 
TDLSDA results for the dipole excitations of the same system. For density modes, the 
numerical calculation reproduces the exact result, and must be considered as a stringent 
test of our three-dimensional (3D) TDLSDA code in view of its application to study other 



excitation modes, like the spin dipole modes considered here, or to address the 
response of more complex systems such as vertical diatomic artificial molecules. 



ar-infrared 
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II. GENERALIZED KOHN'S THEOREM FOR 3D QUANTUM DOTS IN MAG- 
NETIC FIELDS OF ARBITRARY DIRECTION 



In the effective mass approximation, the Hamiltonian for N noninteracting electrons 
confined in a QD by harmonic potentials of freqencies u x ,u y and u z in the x, y and z 
directions, in a constant magnetic field B, is written as 



N 
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where m = m*m e is the electron effective mass, A is the vector potential which we write 

— * — * 

in the symmetric gauge as A — (B A r)/2, and g* is the effective gyromagnetic factor. In 
the numerical calculations we have used the values of to*, dielectric constant e and g* of 
GaAs, namely, to* = 0.067, e = 12.4, and g* = —0.44. Without loss of generality, we write 
B = B(sm9, 0, cos#), where 9 is the angle between the magnetic field and the z axis. Using 
hereafter effective atomic (dot) units (d.u.) h = e 2 /e — to — 1, it can be easily checked that 



i 1 ->2 i 1 / 2 2 , 22, 2 2A , r 

h o = t;P +o \ u x x + u y y +UJ z z ) + hm > 



(3) 



where 



h m = \ool\x 2 cos 2 6 + y 2 + z 2 sin 2 9 - 2xz cos 9 sin 9} + \g* s ^Bf]cT 
-iu c [s\n9[yf z -z^)+cos9[xf y -y£) , 

— * 

and r\ c = +1(— 1) for a =| (J,) with respect of the direction of B. The interacting system is 
described by the Hamiltonian H = H + V, where V is the electron-electron interaction 

N 



V 



(4) 



i<j=l 



The single particle Hamiltonian ho can be exactly diagonalized: 



ho = J2^ a (^c+c a + ^ , 



(5) 



with the creation operator c+ given by 

c+ = a a x + b a y + c a z + i[d a p x + + / Q p z ] . 



(6) 
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The operator c a is the Hermitian conjugate of c+, and the coefficients a a , f a in Eq. © 
are determined by solving the equation 

[h , c+) = u a c+ , (7) 

together with the normalization condition [c~,c+] = 1. One gets the following homogenous 
system of linear equations: 



% 11 
auj + - buj c cos 9 + d(uj 2 + -u; 2 cos 2 9) — -^fu 2 sin 9 cos 9 = 

% % 1 

-aui c cos 9 — bu cu c sin 9 — e(u 2 H — u 2 ) = 

2 2 y 4 c; 

z 1 1 

-6u; c sin — cuj + -rduj 2 sin # cos 9 — f{u) 2 z + -u 2 sin 2 9) = 

a + gLj + — eu c cos = 

b — -duj c cos 9 + eu + - /cj c sin = 

c — -eu) c sin 6 1 + /a; = (8) 
from which the energies wi, 6^2 , ^3 are obtained by solving the secular equation (x = uo 2 ): 



3 2/ 2 1 2 1 2 1 2\ 

I 2 2 2/i 22 2/i 22 22 22 2 2\ 

-x(w c w x COS 9 - uj c uj z COS - UJ y UJ z - U X UJ C - UJ x UJ y - u x u z ) 

-u 2 x u;y z = . (9) 

For each energy solution u> a , Eqs. (jHJ) supplemented with the normalization condition 

Re[ad* + be* + cf*] = -^ , (10) 

where the * indicates complex conjugation, give the coefficients a a ,...,f a . 
Defining = J2f=i c t,a-> ^ * s eas y to prove from Eqs. and (0) that 

[V,C£] = 0, (11) 

which is valid not only for the Coulomb interaction but also for any V that depends only 
on the relative distance between any two electrons. From Eqs. (J7|) and (jllj) it follows 
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immediately that the Hamiltonian H = Y2j ho{j) + V satisfies the equation of motion 



[H, Ct\ = ±u a Ct , (12) 

which implies that if |n) is an eigenstate of H with energy E n so are C^\n) with energies 
E n ±u a . Since in the long- wavelength limit light photoabsorption is induced by the dipole 
transitions, and the dipole operator can be written as a sum of C£ and C~ operators, one 
recovers that dipole transitions can only occur from an eigenstate \n) to the eigenstates 
C^\n), and the absorption spectrum consists of three peaks of frequencies ct>i,o>2, and uj%. 

As an example, we have solved the homogenous system Eqs. (JBJ) in the case of an axially 
symmetric parabolic QD with u x — u y — uj by taking ujq =4.42 meV and u; 2 =18 meV, for 
some values of the 9 angle. The results for the three u a energies are reported in Fig. ^ Due 
to our choice of axial symmetry, at B = the energies lo\,i are degenerate and equal to co>o> 
and co>3 coincides with uj z . At 9 = 0, when B is parallel to the symmetry axis of the dot, the 
solution for o; 1>2 coincides with Peeters' expression Eq. (JTJ), and co> 3 is B independent and 
equal to u z . At 9 = 90°, when B is perpendicular to the symmetry axis of the QD, it is 
possible to write down an analytical expression for u; 2 ,3 : 

"Is = \ + ^ + ± A /(^ 2 + ^ + ^) 2 -4^o 2 } , (13) 
whereas u)\ is B independent and equal to ojq. 
The dipole strength 



n j 

where e is the polarization direction and D^. = ^ . e • fj the dipole operator, can be also 
analytically calculated, and is different from zero only when the excitation energy (E n — Eq) 
is equal to the frequencies u a . One can show that the energy- weighted sum rule mi 



mi 



J S{u)udu = J2( E n ~ E ) \(n \D 6 \ 0)| 2 5(u - (E n - Eq)) = y , (15) 



is exhausted by the three excited states at the energies ^1,2,3- The m_i sum rule, related with 
the static polarizability, can also be worked out yielding a result that is also B independent: 



du N ^ e 2 



S(u) — = — V -|, q = x,y,z (16) 
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(see also Ref. Illl ). By changing the values of 9 and B, some strength moves from one state 
to another, always preserving the values of m_i and mi. 



III. TDLSDA CALCULATIONS 

We now turn to our main goal wich is to numerically study within TDLSDA the dipole 
modes both in the density and spin-density channels. As we have indicated in the Sec. 
I, in this approach the generalized Kohn theorem discussed previously still holds and the 
TDLSDA numerical results provide only a test of the 3D+time TDLSDA code. In the 
spin density channel the theorem does not hold, and the calculations provide an alternative 
prediction for this mode to that provided by the explicit evaluation of the spin-density 



correlation function. 



We recall that spin dipole modes have been experimentally detected 



in Raman scattering experiments. 0, 0] 

To obtain the dipole strength in the two channels, we study the time evolution, following 
an initial perturbation, of the dipole signal 

V(t) = e-(D), (17) 
where D = . r} for density dipole, and D = ^ ■ r}<7| for spin dipole modes, and (D) means 
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average over the time dependent state. This method has been used in Refs. [l^, ll 
to study charge and current modes of two-dimensional QD and quantum molecules; what 
follows is a generalization to the 3D case. 

To calculate T>(t), we firstly solve the static Kohn-Sham (KS) equations 

-\ (& + w + &) + \^ x2 + u yy 2 + 

+V h + V xc + W xc r ]a + h m \ * a (x,y,z) = e a V a (x,y,z) , 
where the expression in the square brackets is the KS Hamiltonian T~Cks- More pre- 
cisely, V H (x, y,z) is the direct Coulomb potential, V xc = 5£ xc (n,m)/5n\ gs , and W xc = 
5S xc (n, m)/5m\ gs are the variation of the exchange-correlation energy density S xc (n, m) writ- 
ten in terms of the electron ground state (gs) density n(x,y, z), and of the local spin mag- 
netization m(x,y,z) = nUx, y, z) — nHx, y, z). The exchange-correlation energy has been 
taken fr „ m pLdew Id zi g e r ,Q J^L) has been exacted a 8 indeed in Re, 



13j . It is worth noticing that if B ^ the single particle wave functions ^/ a (x,y,z) are 



complex and their real and imaginary parts are coupled by h r 



The KS equations are solved in a 3D mesh with A x = A y = 0.51 d.u., and A z = 0.135 
d.u., by using 11 points formulas for the differential operators and a fast-Fourier transform 
to solve the Coulomb potential (for more details see Ref. Q). 

Appropriate static solutions of the KS equations are then used as initial conditions for 
solving the time- dependent KS equations 

i^ <T (*,t)=HKsM*,t), ( 18 ) 

Specifically, to describe the interaction of the system with an external dipole field the gs 
orbitals are slightly perturbed according to 

V£(r) = £/W(r) (19) 

with 



U = exp[zAe-r] , (20) 

for the density dipole modes, and 

U = exp[i\r] a e ■ r\ (21) 

for the spin dipole modes. Eqs. (|T9"j) and (j2*U)l give rise to an initial state in which all the 
electrons of the dot have a rigid velocity in an arbitrary direction e. In the spin dipole case, 
Eqs. (|19j) and (J21|) give initially to spin up and spin down electrons a rigid velocity field 
in opposite directions. The parameter A is taken small enough to keep the response of the 
system in the linear regime. 

We have solved Eqs. (fTHjl following the method of Ref. [3 One writes 

^(r, t + St)- t-5t)= [e' mKS 5t - e m « s 5t ] M*, t) (22) 
and then expands the square bracket in a Taylor series: 

W(r, t + St) = t-St) + 2 ]T (-iy +1 {St H KS f j+1) M^ t) . (23) 

3=0 [ZJ + 

At t = we follow the iteration procedure of Ref. jly up to seventh order. After the first St 
step, we use Eq. (J2HJ), taking in the expansion j max = 3. Typical St are ~2x 10~ 2 d.u. (1 



d.u. is ~ 5.55 x 10~ 14 s). After 15000 iterations the total energy is conserved with a relative 
error smaller than 10~ n . The dipole signal Eq. (j!7|) is calculated every time-step. 
In the linear regime the Fourier transform of the dipole signal 

V(uj) = J dte lwt V(t) (24) 
is directly related to the dipole strength S§(a>) Eq. (fM"|) by 



S 6 {u) = \V(cu)\ . (25) 

Hence, a frequency analysis of T>(t) provides the absorption energies and their associate 
intensities. jl7] 

A real time simulation of the dipole evolution is shown in Fig. El for a, N = 6 quantum 
dot for 9 = 22.5° and B — 5 T. The analysis is made for the three components of the 
dipole operator in the density channel, and the frequency analysis of the numerical signals 
providing the density dipole strength is plotted in Fig. El for different values of 9 and B. In 
this figure, the energies uj\ ;2 ,3 yielded by the exact calculation of Sect. II are indicated with 
arrows. The peak energies of the TDLSDA strength are also reported with circles in Fig. 
[0 From these figures one may conclude that our 3D TDLSDA calculations reproduces very 
accurately the exact results of the previous section. From Fig. El one may also see how, by 
changing the values of 9 and B, some strength moves from one mode to another. 

Performing the same simulation in the spin channel we obtain the results shown in Figs. 
E] and Differently from the density channel, in the spin channel all the strength is mainly 
concentrated in one collective low energy mode. The peak energy is lower in the spin than 
in the density channel. This is due to the character of the TDLSDA particle-hole residual 
interaction, which is attractive but weak in the spin channel, and repulsive and rather 
strong in the density channel, shifting the TDLSDA response from the single-particle one in 
opposite directions. [6] Increasing the value of B, the spin mode becomes progressively softer. 



o 


: a 


> 


22 



221 ] . A spin mode with this feature has been observed 



in GaAs quantum dots|4| using Raman spectroscopy, as well as in quantum wells 23] in a 
perpendicular magnetic field. In particular, in the quantum well experiment, evidence for 
the spin instability at some values of electron density and B has been reported. [23] The 
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agreement between our calculation and the experimental findings of Schiiller et al. y in 
QD's is qualitatively good for low B values (see also Ref. 0). The lack of experimental 
observation of collective spin states at larger B values might be considered the signature of 
the above mentioned instabilities. However, it cannot be excluded that Landau damping, 
which in this energy range is particularly strong, prevents its experimental observation. 
Besides, correlations beyond TDLSDA which are important for non-integer filling factors, 
might also quench the spin mode. 



IV. SUMMARY 



The dipole and spin dipole responses of a 3D quantum dot in a magnetic field of arbitray 
direction have been analyzed by means of real time simulations within TDLSDA. The frag- 
metation of the strength due to the non circular shape of the QD and to the 9 angle between 
the magnetic field B and the axis perpendicular to the plane of motion of the electrons has 
been studied. The density dipole strengths splits in three collective states which are exactly 
found at the energies predicted by the generalized Kohn theorem. By changing the values of 
9 and of B, some strength moves from one state to another in a way which could be easily 
detected. 

The spin dipole strength is mainly exhausted by a single, soft collective mode whose 
energy goes to zero as B increases, and as the direction of the applied magnetic field goes from 
perpendicular to parallel to the plane containing the dot. For some values of B, TDLSDA 
predicts spin instabilities similar to these observed in Raman spectroscopy experiments. 

A key result from our study is the circumstantial evidence that it is possible to develop 
a 3D code implementing TDLSDA with a magnetic field in an arbitrary direction with very 
high accuracy. This opens the possibility of addressing the far-infrared response of more 
complicated systems, such as double QD's[ll3] or quantum rings vertically coupled, whose 
description using the density-density or spin-density correlation functions [6( is prohibitive. 
Work along this line is now in progress. 
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FIG. 1: Exact dispersion relation for dipole density modes. Filled circles are from TDLSDA 
calculations 

FIG. 2: Real time evolution of the dipole density signal. Shown are the x, y and z-components of 
the dipole signal for 6 = 22.5° and £=5T. 

FIG. 4: Real time evolution of the dipole spin signal. Shown are the x, y and z-components of the 
dipole spin signal for 9 = 67.5° and 5= IT. 



FIG. 5: Strenght function (arbitrary units) for the dipole spin response corresponding to different 
angles and magnetic fields. 



FIG. 3: Strenght function (arbitrary units) for the dipole density response correponding to different 
angles and magnetic fields. 
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